A new model for magnetoreception 
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Certain migratory birds can sense the earth's magnetic field. The nature of this process is not 
yet properly understood. Here we offer a simple explanation according to which birds literally 
'see' the local magnetic field: Our model relates the well-established radical pair hypothesis to the 
phenomenon of Haidinger's brush, a capacity to see the polarisation of light. This new picture 
explains recent surprising experimental data indicating long lifetimes for the radical pair. Moreover 
there is a clear evolutionary path toward this field sensing mechanism: it is an enhancement of a 
weak effect that may be present in many species. 



It is well established that certain migratory birds can 
detect the direction of the Earth's magnetic field, and 
use this as a navigational compass.^ An obvious expla- 
nation for this remarkable ability would be the use mag- 
netised materials in the bird's body, which change orien- 
tation with the external field. However for species such 
as the European Robin the evidence points to a very dif- 
ferent mechanism: the prevailing hypothesis is that field 
orientation is initially detected through its influence on 
photo-excited electronic spins. Photons are evidently im- 
portant since the birds can only orientate in a magnetic 
field when it is light, ^ ^ and with an undamaged visual 
system.!^ Meanwhile a recent observation provides strong 
support for the role of electron spins: the birds' are dis- 
orientated by a weak oscillatory field whose frequency is 
close to the resonant frequency for an electron in the 
Earth's magnetic field. ^ ^ While these detailed results 
come from bird studies, there is evidence of a similar sen- 
sitivity in the American cockroach, suggesting that this 
form of magnetoreception may occur in diverse species.— 

These findings have led to the popularity of the Rad- 
ical Pair (RP) model,'^'^^!^ which begins to explain how 
light activated magnetic sensing could happen. The cen- 
tral feature of the model is that the optical excitation 
of certain biomolecules leads to there being a fraction 
of molecules in a spin triplet state, with a spatially- 
separated pair of spins. In magnetically anisotropic sys- 
tems, the number of spin triplets depends on the orien- 
tation of the magnetic field. If the bird can detect the 
population of spin triplets somehow, then an optically 
activated avian compass is possible. 

The route by which electron spin states translate to a 
macroscopic signal is inadequately understood. A typical 
explanation is that some signature chemical is synthe- 
sised only when the triplet state decays. This chemical 
may then interfere with the normal process of vision, or it 
might be detected by some independent sensor structure 
in the eye. Explanations of this kind have two short- 
comings. They involve a complex chemistry which must 
somehow have evolved within the eye, independent of 
(but consistent with) the process of normal vision. More- 
over, this model would seem to function best when the 



cycle time, i.e. the time for production of the signature 
chemical(s), is short - shorter cycles would lead to higher 
rate of production and thus better signal-to-noise ratios. 
However, in the real system it seems that the opposite 
is tr ue: t he relaxation time of the RP is extraordinarily 

long.isini 

Here we will describe a model of the compass in which 
it is straightforward to understand the need for long 
lifetimes; indeed this property is so crucial that the 
molecules involved could have evolved through natural 
selection of slow electronic decay rates. We will further 
describe why it is possible to do away with any apparatus 
to detect chemical products altogether. In essence, the 
Earth's magnetic field translates to a local electrostatic 
(or strain) field which directly modulates vision. This 
may be seen as an evolutionary enhancement of an inher- 
ent sensitivity, analogous to the well studied 'Haidinger's 
brush'. We believe that our new version of the RP model 
is thus a simpler and more complete hypothesis, while 
being more consistent with experimental findings. 

Our model borrows the successful features of the RP 
model, also requiring light to create charge separated 
radical pair triplets. It is natural to assume that such 
optically created triplets are formed somewhere in the 
eye, either in lens system or in the retina. The retina 
seems more likely, as it is already integrated with a sys- 
tem to initiate signals to the brain in normal vision. In 
our mechanism, radical formation would be indeed be 
in the retina, probably in the most sensitive region (the 
macula lutea) directly opposite the lens. The magnetore- 
ception system then 'hitches a ride' on the normal vision 
system, modulating it by the electric dipole of, and 
therefore electric field created by, the triplet state. 

We base our model on a rod-like molecule with the 
following three properties. First, it possesses an ex- 
cited state singlet to triplet branching ratio dependent 
on an external magnetic field. This phenomenon has 
been observed in, for example, self-trapped excitons in 
alkali halides^^ and NV~ centres in diamond. Second, 
we require a charge-separated excited triplet state with 
a sufficient lifetime to influence the visual process in the 
retina. Sufficiently long triplet lifetimes of a few millisec- 
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onds are not uncommon. Third, we require a spin level 
structure that allows the long-lived triplet state to be 
dephased by a resonant RF field. We shall discuss one 
possibility for this later. Further, in Fig. [T] we display 
how these molecules might be arranged in a circular pat- 
tern on the surface of the retina. This kind of pattern has 
already been proposed for lutein molecules as a possible 
explanation of the fact that some people can directly see 
light p olariz ation (a phenomenon known as Haidinger's 
brush) JEli6J ^j^g alignment could come from the known 
radial orientation of nerve fibres. The molecules in our 
model would be magnetically anisotropic, not unexpected 
for a rod-shaped molecule. 

Optical excitation gives rise to the formation of charge- 
separated triplets states for certain orientations of the 
magnetic field with respect to the molecular axes - and 
the result is an electric field distribution on the retina 
which refiects the orientation of the magnetic field. Rel- 
atively weak electric fields between 10^ — 10^ V/m affect 
the photoenergetic reactions of bacteriorhodopsii^^^^^ as 
well as the cis-to- trans isomerization of many other com- 
plex molecules .I^^M^ An electron-hole dipole with average 
charge separation of one nanometre produces an electric 
dipole field with magnitude 10^ V/m up to a distance of 
10 nm, while a field of the order of 10^ V/m even extends 
to 25 nm. Each compass dipole thus possesses a sizeable 
'sphere of influence' around it in which it could directly 
affect the photoisomerization of retinal, which forms part 
of rhodopsin^^ - and the bird would literally be able to 
'see' the magnetic field as a superimposed feature in its 
normal visual image. 

Let us now consider an example system that would 
exhibit the features required by our model. We will de- 
scribe the simplest possible molecular energy level struc- 
ture that would work, though it is of course likely that 
any real system will have extra features. The scheme 
we have in mind has four important energy levels and is 
sketched in Fig. [2] 

We imagine that our biomolecule, like most others, has 
a singlet ground state l^o). Light can then excite pop- 
ulation to a higher lying singlet \Si), probably through 
a cascade of non-radiative decays of other higher excited 
singlets. Such transitions would be strongly allowed by 
dipole selection rules. The excited singlet would nor- 
mally simply decay back to the ground state, but it is 
also possible that population could go into a (degenerate 
or close to degenerate) triplet state |Ti). As we discuss 
below, the rate of this interconversion - or intersystem 
crossing (ISC) rate - can be dependent on the geomag- 
netic field, and lies at the heart of the magnetoreception 
mechanism. Finally, |To), reached through optical de- 
cay of |Ti), is a lower- lying, long-lived triplet state with 
a charge- separated orbital configuration of electron and 
hole. Any population in this level has an associated elec- 
tric dipole moment, and this triggers the visual stimulus 
for the compass. 

The dependence of the intersystem crossing rate on the 
geomagnetic field must be associated with an anisotropic 




excitation cycle of the step in the vision 
compass molecule process 




FIG. 1: a: Left: Schematic diagram showing how molecules 
are aligned on macula lutea; a combination of light and a mag- 
netic field could induce dipole moments for certain molecular 
orientations. Right: These dipoles would create an electric 
field which would allow the bird literally to 'see' the magnetic 
field direction, b: Cycle of the compass molecule: following 
photoexcitation from |5'o) to \Si), the branching ratio of di- 
rect relaxation into the ground state or via a long-lived triplet 
state I To) depends on the orientation of the molecule with the 
geomagnetic field. The purple colour of |To) denotes a charge 
separated state with an electric dipole moment, thus affect- 
ing the isomerization of retinal, which is a crucial step of the 
visual process. 




FIG. 2: Schematic diagram showing the key features of the 
compass mechanism: the system relaxes back into the singlet 
ground state \ So) after most photoexcitation events. However, 
there is a small intersystem crossing rate, which depends on 
the orientation of an asymmetric ^-tensor with the geomag- 
netic field. Population in |Ti) relaxes into a long-lived triplet 
state I To). Transitions between different charge states are de- 
picted in red, while a change of the spin state is shown in 
blue. The decay from |To) to |5'o) involves both a charge and 
a spin transition. 
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term in the Hamiltonian. The origin of this term is not 
important; a hyperfine couphng between the electron in 
the opticahy excited exciton and a nu clear spin has been 
widely proposed in the literature.^ In order to keep our 
discussion as simple as possible, we will assume that the 
excited triplet state |Ti) has an isotropic electron ^-factor 
= 2 and a uniaxially anisotropic hole tensor, 

(2^5g \ 

^h=\ 2^5g , (1) 
\ 2 ^Ag J 

and as specific examples, we assume Sg = and Ag = 
0.2. We have checked that the qualitative predictions 
of our model also work for the case of an anisotropic 
hyperfine coupling similar to the one described in Ref. JTI 
Following the photoexcitation, the system relaxes to the 
pure singlet state \Si) which is degenerate with the triplet 
level I To). We write for the Hamiltonian at this stage of 
the process: 

Hisc = \^JiB {9eB • Si + B • • S2) , (2) 

where /i^ is the Bohr magneton, B is the magnetic field 
vector, and = (cr^^, cr^, cr^)* is the spin operator for 
electron (i = 1) and hole (i = 2). The factor 1/2 
accounts for the fact that all our Pauli matrices have 
eigenvalues ±1. The magnetic field strength in Frank- 
furt (the site where the relevant experiments where per- 
formed"^^) is ^0 = 47 /iT.I^The field's orientation with 
respect to the ^-tensor is determined by the angles 6 and 
0, B = 5o (cos (j) sin ^, sin sin ^, cos 0) . 

Based on Hamiltonian we obtain the following ma- 
trix elements for the three triplet sublevels = | tt)^ 
l^o) = (I n) + I m/V2 and = I U) of |Ti): 

{Si\Hjsc\to) = fiBBoAgcose, (3) 
{Si\Hisc\t±) = ±fiB Bo Sg sin ee^"^/V2. (4) 

To obtain a signal that depends on the relative orienta- 
tion of ^-tensor and field, as is required for a compass, we 
must thus have 5g ^ A^, a condition which is fulfilled by 
our particular choice of parameters. Owing to the axial 
symmetry of Hamiltonian ([2]), ^ is unimportant and we 
need only consider G [0,7r72]. 

We now use a phenomenological Lindblad master equa- 
tioiP^ to model the evolution of the density matrix that 
describes the quantum dynamics of our (open) molec- 
ular system. The optical excitation between l^o) and 
\Si) is modelled as an incoherent process with a Lind- 
blad operator Px = \Si){So\ with associated rate jx- 
Similarly, the decay events are described by Lindblad op- 
erators Ps = \So){SiIPt = |To)(Ti| and Pq = |^o)(To| 
with respective rates 75, 7^ and 70 as depicted in Fig.|2j 
Using only the matrix element Eq. (|3| as the effective 
Hamiltonian H and all of the above Lindblad operators, 
we obtain as the master equation governing the time evo- 




FIG. 3: Steady state population of |To), T, is shown as a 
function of the state's lifetime I/70 and the angle between 
the axial ^^-tensor and the Earth's magnetic field for = 0.2. 
For other parameters see main text. The 2D-plot in the upper 
right follows the cos^ dependence of the squared relevant 
matrix element [see Eq. (|3|]. 



lution of the system's density matrix p(t) p5 | 26 ] 
p=-'-[H,p] + Y,li (PipPi - I (p}PiP + pPlP^) ) -(5) 

We are interested in the steady state population T of 
the charge- separated triplet level |To), which is found by 
setting the LHS of Eq. ([5| to zero and is given by 

^lx9{0) 

^1X9{0) + 7or ((75 + 7x)^2 + 4 (1 + ^) 9{0)) 

(6) 

where F = (75 + 7t) and g{e) = \{Si\Hisc\to)f ^ In the 
regime of interest, the lifetime of |To) is much longer than 
that of the other excited states, i.e. 70 <C 75 ~ 7t. It 
is also reasonable to assume that jx <C 75 and 75 
jiBB^Ag /h^ and to a good approximation we therefore 
find T oc g{0)/jo. 

Fig. [3] shows a full numerical solution of the triplet 
population T as a function of the angle 6 and the triplet 
lifetime I/70. Here, we have assumed the lifetimes of 
the \Si) and |Ti) states are 1 ns, and used an excitation 
rate 7^ = 10^ s~^. As shown in Fig. [3]a longer excited 
lifetime would increase the average electric field (i.e. the 
product of the molecule's dipole moment and the num- 
ber of dipoles present at any time), but importantly, it 
also gives each individual dipole more time to have an 
effect on other processes in the vision system such as the 
isomerization of any nearby retinal. 

Recent experiments show that a very small oscillating 
magnetic field can disrupt the bird's ability to orient.^El 
In Ref. 6, the authors report that a perturbing magnetic 
field of frequency of 1.316 MHz (i.e. the resonance fre- 
quency of an electron spin with a ^-factor of 2) with a 
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field strength of only 150 nT suffices to completely dis- 
orient the birds. The authors note that the MHz fre- 
quency immediately implies a bound on the time of the 
process (since the field would appear static for a suffi- 
ciently rapid decay back to the molecule's ground state 
in less than a microsecond). Moreover, considering the 
oscillating magnetic field strength implies a much longer 
process time of at least several hundred microseconds to 
give the weak radiofrequency field sufficient time to affect 
the spin state. 

There are many possible mechanisms by which a weak 
RF field could disrupt the compass mechanism by short- 
ening the lifetime of such a long-lived triplet state. In the 
following, we discuss one simple possibility. We focus on 
the lower-lying triplet to ground state transition. For un- 
derstanding the directional sensitivity we needed no more 
than those two levels, but in order to understand the res- 
onant effects of a magnetic field we must now explicitly 
include the distinct spin states, see Fig[4j We distinguish 
the specific triplet states |to), and we label the 
corresponding singlet level as \S^). This latter level will 
have a fast, spin-allowed decay to the molecule's ground 
state |5'o). State l^') thus separates the spin and the 
charge transition of this relaxation process. Importantly, 
the auxiliary level could be eliminated from the dynam- 
ics so long as the decay 75/ is large enough, reducing 
the model once more to the simpler picture displayed in 
Fig.H 

In our illustrative example, we assume an electron g- 
factor of = 2 and = 2.2 for the hole (both isotropic). 
Working in the basis where the z-axis is defined by the 
applied static magnetic field direction, the triplet now has 
the three sublevels shown in Fig. |4j and the Hamiltonian 
in an RF field reads: 

^RF = ^I^B [Bo {geCTzl + ghCrz2) 

-^Brf COS Ujt {geCTxl ^ 9hcrz2) ]. (7) 

Here, the oscillatory field of strength 5rf = 150 nT and 
frequency uj is applied orthogonal to the static field 
since only its perpendicular component affects the com- 
pass.f^ Note that the three sublevels |to) and 

of Fig. [4] take a different form to the sublevels of |Ti) 
that we have considered before, whose axis was defined 
by the anisotropic ^-tensor of the hole. (Following a spin- 
preserving relaxation from |Ti) to |To), the state within 
T is cosi9|to) +sin(9(|t+) - \t_))V2.) 

Fig. [5] shows the surviving |To) population as a func- 
tion of time, obtained with a master equation that is 
isomorphic to Eq. ([s]), but now based on Hamiltonian 
with the Lindblad operators P5/ = |5'o)(5"| for the 
decay, and Px^i = {crx)i for the relaxation of electron and 
hole spin, respectively. The decay rate is 75/ = 10^ s~^ 
and we assume slower but equal spin relaxation rates of 
0.2 X 10^ s~^. As the initial state, we assign half of the 
population to \to) and the other half is equally distributed 
between the sublevels and consistent with the 
state resulting from a |Ti) decay. 



^ spin relaxation ^ ^ 




FIG. 4: Possible explanation of the relaxation from |To) to 
the ground state \So): the charge relaxation follows a spin 
transition to the auxiliary singlet level \S'). Population in 
I to) is connected to \S') via an ISC and thus decays quickly, 
whereas population in \t±) is trapped for the duration of the 
spin coherence time. However, a resonant RF field mixes the 
spin states, leading to faster relaxation of the entire triplet 
population. 
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FIG. 5: Surviving |To) population as a function of time and 
RF frequency uj. After an initial fast decay of the |to) popula- 
tion (see text), a slowly decaying plateau of triplet population 
is reached. However, on resonance with the either the elec- 
tron or the hole spin, the oscillatory field drastically shortens 
the triplet lifetime. See text for parameters. 



The pronounced kink at short times in the data of 
Fig. [5] shows how the |To) relaxation proceeds in two 
stages: population in \to) undergoes a direct ISC to l^*'), 
subsequently decaying to the ground state in much less 
than a millisecond. The other half of the triplet popu- 
lation in the \t±) sublevels survives much longer, until 
all spin states are eventually mixed by the Px^i = {crx)i 
spin decoherence processes. Only when the oscillatory 
field is resonant with either the electron or hole spin (at 
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1.316 or 1.447 MHz due to the different ^f-factors), are 
all triplet levels mixed on a sub- millisecond timescale, 
leading to a fast decay of the entire triplet population. 
The corresponding two valleys in the data elucidate how 
the resonant RF field thus severely reduces the excited 
triplet lifetime. Such a reduction of the |To) lifetime 
would plausibly affect the compass mechanism and dis- 
orient the bird. 

We note that the general behaviour displayed in Fig. |5] 
does not depend on the specific choice of parameters, so 
long as the hierarchy of the various processes is preserved. 
Indeed, one can also think of entirely different physical 
mechanisms for the |To) relaxation that would be equally 
consistent with experimental observations. Crucially, all 
possible explanations depend on an excited state with 
a lifetime of more than a millisecond, as the oscillatory 
field strength is simply too weak to significantly affect 
the spin state in a shorter duration. 

The model we have proposed could be tested, and we 
now discuss some possible experiments that might be per- 
formed. 

First, it would be very interesting to probe the RF 
disruption mechanism further. In particular, if one could 
test the ability of a bird to navigate in an RF field across 
a range of static fields, it should be possible to obtain the 
width of the triplet resonance line that causes the disrup- 
tion. Indeed, we would expect to find two resonances for 
the specific model we have proposed - one corresponding 
to electron resonance and the other to hole resonance. 

A second possibility is that birds using the mechanism 



we propose would be very sensitive to the polarization of 
the light used to induce the magnetoreception. The rod 
like molecules discussed in this paper would be expected 
to have a highly anisotropic electric susceptibility, and so 
the method we propose might only work in polarized light 
for molecules of certain orientations. It may be that the 
birds can adapt to this, but by changing the polarization 
of the light to which the birds are exposed periodically 
any adaption could be prevented, and again we might 
expect to see disruption of magnetoreception. 

We have proposed a comprehensive model that would 
allow a bird to sense the direction of the Earth's mag- 
netic field. It relies only on processes that are common 
to standard ideas of vision mechanism and processes that 
exploit perfectly standard features of small biomolecules; 
nothing exotic is required. Further, we have shown that 
a long lived triplet, essential for understanding the ob- 
served disruption of the effect by very weak radio waves, 
could have evolved through natural selection. Interest- 
ingly a number of human biomedical disorders are at- 
tributed - with many doubts and reservations - to low 
intensity, oscillating, electromagnetic fields in the same 
freque ncy ra nge that disrupts magnetism-based bird nav- 
igation.^^^ 
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